Localized Orbital Description of Electronic Structure 
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We present a simple and general method for construction of localized orbitals to describe electronic 
structure of extended periodic metals and insulators as well as confined systems. Spatial decay of 
these orbitals is found to exhibit exponential behavior for insulators and power law for metals. 
While these orbitals provide a clear description of bonding, they can be also used to determine 
polarization of insulators. Within density functional theory, we illustrate applications of this method 
to crystalline Aluminium, Copper, Silicon, PbTiC>3 and molecules such as ethane and diborane. 

PACS numbers: 71.15.-m 71.20.-b 71.23.An 



Localized orbitals, such as Wannier functions^], have 
played an important role in many aspects of quantum me- 
chanical description of electrons in solids and molecules. 
They highlight atomic character of electrons and elu- 
cidate the nature of bonding. Their spatial localization 
makes them useful as basis functions for efficient cal- 
culations of electronic structure that scale linearly with 
system size0, 0, They are often useful in the con- 
struction of lattice model Hamiltonians for structural 
transitions Q as well as strongly correlated systems^- 
Wannier functions are intimately related to geometric 
phases of Bloch electron^, [t| and electric polarization 
of insulators|0,^|. Recent surge of interest in Wannier 
functions is due to their use in modeling transport [f^ | 
and electric field dependent properties of periodic solids. 

Many schemes for construction of localized or Wannier 
orbitals have been presented over the last four decades. 
To name a few, early works of Foster and Boys^| were 
based on maximizing the dipole moment matrix ele- 
ments between orthonormal orbitals making their cen- 
troids maximally apart from each other, followed by 
Reudenberg's |14j approach of maximizing the orbital 
self-repulsion energies. Kohn presented a schemeQ for 
arriving at exponentially localized real orthonormal WFs 
through a variational principle for total energy. Maxi- 
mally localized Wannier functions (MLWF)|f5| were in- 
troduced by Marzari and Vanderbilt that minimize the 
variance of position operator r, in a given subspace of 
states. For example, orbitals maximally localized in 
a spatial direction a (hermaphrodite orbitals ^(j) are 
eigenf unctions of position operator r a in the given sub- 
space. As different components of f do not commute 
when projected into a subspace, MLWFs are obtained 
through numerical minimization of variance of r used as 
a measure of localization. MLWFs have been generalized 
to the cases of entangled bands and to atom-centered 
orbitalsflij. 

WFs are Fourier transform of Bloch functions: W(r) — 
J dke lk r ^k(r). Due to freedom in k-dependent phase 
factor accompanying ifa, WFs are non- unique. Smooth- 
ness of these phases determine localization properties of 



WFs. MLWFs in 1-dimension can be obtained simply 
(with out a variational calculation) through construction 
of Bloch states that are smooth (generated using parallel 
transport along k) and periodic in fc— space 9J. General- 
ization of this idea to three dimensions is not readily pos- 
sible, as the smoothness and periodicity of Bloch states 
as a function of fc Q achieved through parallel transport 
along one direction is disturbed by that along another 
direction in fc— space. Here, we provide a solution by 
connecting Bloch states at different k's by parallel trans- 
port along paths that run outside the k-space. Result- 
ing Bloch functions are periodic and optimally smooth in 
all directions in fc— space, whose Fourier transform yields 
well localized Wannier functions. It is efficient and sim- 
ple because it avoids a variational calculation and issues 
of local minima. Its versatility is demonstrated through 
applications to insulators, metals and molecules. 

In the first step of this method, we develop an aux- 
iliary subspace of highly localized orthonormal orbitals 
with desired symmetry properties described by (a) center 
of the orbital (Wyckoff site) {r K } and (b) irreducible rep- 
resentation (irrep) of its site symmetry group according 
to which it transforms 0, 0] • This choice can be guid ed 
through symmetry analysis of Bloch states 0, lil l2o| in 
the physical subspace of occupied states. For simplicity, 
we use Gaussian form for the radial part and a spherical 
harmonic corresponding to the irrep of localized orbitals: 

*„(R,r) = (a/27r) 3 / 2 e -«l r — ^Y lm (r-^-R) (1) 

where a determines the width of the Gaussian, R is a 
direct space lattice vector and /x is an orbital index, a 
is chosen to be large enough to keep orbitals in neigh- 
bouring unit cells orthogonal. Cell periodic part of 
Bloch functions spanning the auxiliary subspace is given 
by Fourier transform: 
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We note that { (r|i; M k) } are smooth as a function of 
k and satisfy |w M k) = |w M k+G >, G being a reciprocal 
lattice vector. 



2 



In the second step, we perform a unitary transforma- 
tion of Bloch states in the physical subspace (eg. en- 
ergy eigenfunctions of occupied states) such that the open 
path[2l| non-abelian geometric phases between states in 
physical and auxiliary subspaces at fixed k vanish. Over- 
lap matrix S* k n = (w M k|w„k), where \u n k) is cell periodic 
part of energy an eigen state, relates to geometric phase 
matrix Tk through: 

S k = i? k e jI \ (3) 

where R is a positive definite Hermitian matrix. Determi- 
nation of R and I\ is accomplished using singular value 
decomposition of S: 

5 k = tf k £ k V k t (4) 

where Sk is a diagonal matrix with non-negative singular 
values in its diagonal. Vanishing of singular value(s) at 
some k signals a non-optimal choice of symmetries of 
localized orbitals in the auxiliary subspace and the need 
for correction. Using EqnQ] -Rk = U^kU^ and e lTk = 
[/kUj. A unitary transformation on the {u ra k} given by 
a matrix Mk = (U-^V^)* transforms the overlap matrix 
5 k to a hermitian form, 

corresponding to vanishing geometric phases. 

A sketch in terms of paths that connect physical and 
auxiliary subspaces (shown in Fig^) enables understand- 
ing of our method. Closed paths parametrized by k, such 
as C, in the physical subspace can have nontrivial geo- 
metric phases which forbid construction of smooth func- 
tions by parallel transport along k within the physical 
subspace. In the present method, we make use of paths, 
labeled as P, to connect states at two k's in the physi- 
cal space. For k = constant segments of path P (dashed 
lines), geometric phases are made to vanish. Secondly, 
geometric phases along any closed path parametrized by 
k in the auxiliary subspace vanish. Hence, Bloch func- 
tions in the physical subspace transformed with Mk have 
a single valued and periodic phase as a function of k. 

We include occupation numbers while transform- 
ing energy eigenfunctions and obtain 

|u M k) = E M tij\uj*)f„k, (5) 

i 

which are periodic in k by construction and the random 
phases accompanying {it n k} are filtered out in this pro- 
cedure. Localized orbitals are obtained by Fourier trans- 
forming {w pk (r)}: 

< r|* M (R)> = 7^)3 E e * (k ~ R) ' r < r '^ k > ( 6 ) 

where fi is the unit cell volume. Localization of |$ M (R)) 
is determined by the "smoothness" of {u M k(r)}, which 
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FIG. 1: Two horizontal planes are used to display physical 
and auxiliary subspaces, in which the member functions are 
parametrized with k x and k y . Dashed lines indicate paths 
connecting states at a given k in physical subspace with those 
in the auxiliary subspace. 

depends on the precise location of centre and symmetry 
properties of localized orbitals in the auxiliary subspace. 
A quantitative idea of the smoothness is obtained with 
the Berry connection matrix B a (k): 

B^(k) = -7m<«^|^-|tU>, (7) 

and its Fourier components. If the matrix B a {k) is diag- 
onal and its diagonal entries are constant as a function 
of k a (it has vanishing Fourier components for R/ 0), 
$ M are maximally localized in a direction^ 0]. For the 
systems studied in this work, diagonal elements of B are 
found to be constant as a function of k, indicating good 
localization properties of 

We note that functions |$ M (R)) are not maximally 
localized 15] by construction. However, if desired, the 
MLWFs can be readily obtained through a single step 
of maximal joint diagonalizationj^ of the three compo- 
nents of position operators in <3>-basis: 

(ro.W.Jii.ife = / **(Ri,r)r Q $^(R2,r)dr, (8) 

J(lN ka 

Electronic part of electric polarization is obtained using: 

Pa = £ E<M R = °)M<MR = o)). (9) 

Inclusion of occupation numbers in Eqn [S] ensures a 
simple and general form for density matrix p(r, r') = 
J2r /i ( r l*M(-ft))(^V(-ft)l r ')- As the present analysis at 
k is decoupled from that at other k's, it applies equally 
well to simple, compound and entangled bands, and can 
be performed efficiently on a parallel computer. For 
molecules or confined systems treated within periodic 
boundary conditions, it can be readily applied using a 
single k = (000) point; it equally well applies to con- 
fined systems treated with open boundary conditions. 
In the pre sent scheme, construction of atom-centered 
orbitals [lq for covalent systems necessiates expansion of 



3 




~^T/ — ~ 

► 50% 













> 




■ ( 



FIG. 2: (a) and (b) shows 2p y and 2^ orbitals respectively of 
a bridging oxygen in PbTiOs. (c) shows the nature of decay 
of the 2p x orbital shown in (b). (d) shows a bond centered 
orbital representing the Si-Si a-bond in bulk Si. 



the physical subspacc to include anti-bonding bands. It 
is somewhat similar to physically intuitive derivation of 
Wannier functions in Ref 23. 

We illustrate our method within density functional 
theory through applications to insulators, metals and 
molecules. The energy eigen states, which are the main 
inputs to our method, are calculated using ABINIT 
implementation |24| of density functional theory and 
norm-conserving pseudopotentials. We use Monkhorst 
Pack meshes (finer than a 9x9x9 mesh) of k-points to 
sample Brillouin zones. Isosurfaces of localized orbitals 
(indicating isovalue as percentage of its maximum in fig- 
ures) have been generated using a visualization software 
XcrysDenjUJ. 

For PbTiC>3 in the cubic perovskite structure, symme- 
try of occupied bands implies the choice of atom-centered 
orbitals: Pb-centered (s— and d— like), Ti-centered (s— 
and p— like) and O-centered (s— and p— like). Oxygen 
centered localized orbitals with p symmetry form two 
groups: (a) ones perpendicular to the -O-Ti-O- chain 
(FigEJ(a)), (b) ones along the -O-Ti-O- chain (FigEI(b)). 
A former has it— like overlap with d xy orbital of Ti and 
a variance of 1.44 A 2 , while the latter has a— like over- 
lap with d x 2 orbital of Ti and a smaller variance of 1.14 
A 2 . For an O-centered orbital with s symmetry, we 
get a variance of 0.54 A 2 (compared with 0.52 A 2 of 
Ref. 26]). We determined Born effective charge (Z*) 
of Ti using Z* = ilAP a / Ad a , where Ad a is a small dis- 
placement of Ti atom that changes the net polarization 
(P* = P l a on + Pa) by AP a . Our estimate of Z*(Ti] 
7.02 agrees well with a linear response calculation)! 
with contribution from each O-centered p orbital in group 
(a) and (b) of 1.6 and 0.8 respectively. Our estimate of 



Z*(Pb)=3.86, also in good agreement with the linear re- 
sponse result, has a contribution of 1.4 from O-centered 
orbitals in Ti-0 planes. Spatial decay of an O-centered 
orbital in group (b), shown in Fig. Etc), exhibits a power 
law times times exponential form|2|. The power-law ex- 
ponent obtained from a fit is —0.75, consistent with a 
theoretical prediction!^. We confirmed that the diag- 
onal elements of the Berry connection matrix (Eqn 7 do 
not change with k, indicating individually maximal lo- 
calization of these orbitals. 

For Si in diamond structure, a covalent semiconduc- 
tor, we illustrate construction of bond and atom cen- 
tered orbitals corresponding to different choices of sub- 
spaces. Symmetries of Bloch states in the occupied sub- 
space (D=4) dictate a choice of bond-centered orbital 
with full site symmetry. Corresponding bond centered 
orbitals (shown in FigEJ(d)) have a variance of 2.2 A 2 
(compared with 2.05 A 2 of an MLWFpjJ) each. We find 
that the diagonal elements of Berry connection matrix 
are indeed constant as a function of k 7 implying maximal 
localization of an orbital individually. In the construc- 
tion of atom-centered orbitals we use double the num- 
ber of bands in physical subspace, and atom-centered 
orbitals with s and p symmetry. Our scheme provides 
two options: (a) treat silicon as a metal (using occupa- 
tion numbers in EqnEJl, which generates nonorthonormal 
atom-centered orbitals which exactly span the occupied 
subspace, and (b) not use occupation numbers, which 
generates orthonormal atom-centered orbitals that span 
the occupied as well as unoccupied states. The former 
can be important in studies of bonding, whereas the lat- 
ter could be useful as basis in O(N) methodology. 

We illustrate localized orbitals in metallic systems 
characterized by multi-centered bonding, ie. sharing of 
electrons among more than two atoms, with examples of 
Cu and Al in the FCC crystal structure. Due to partial 




FIG. 3: Localized orbitals of metals: an orbital centered at 
the tetrahedral site of Al (a), and its spatial decay (b); atom 
centered LO with 3d z 2 symmetry (c) and an LO centered at 
the octahedral site of Cu (d) . 



4 




FIG. 4: (a) shows a bond centered orbital representing the 
C-C u-bond in C2H6. (b) shows the B-H-B three center bond 
in B 2 H 6 . 



to extended metals and insulators as well as finite sys- 
tems. While demonstrated here for electronic problems, 
it can be readily used in treatment of phonons. It should 
be useful for a large number of problems in condensed 
matter physics. 
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occupancy of bands in metals, the dimensionality of aux- 
iliary subspace is taken to be greater than largest number 
of occupied bands at any k, the same as the number of 
energy bands calculated in DFT calculation with temper- 
ature smearing used for occupation numbers. In the case 
of Al, we used (guided by symmetry and singular values) 
spherically symmetric orbitals centering at the tetrahe- 
dral sites (1,1,1)^ in addition to atom-centered s— like 
orbital to construct an auxiliary subspace. From nor- 
malization, we find each localized orbital centered at the 
tetrahedral site (Fig|3J(a)) and the atomic site to be occu- 
pied with 1.2e and 0.6e respectively. The former provides 
4— centered directional bonding among the four Al atoms 
equidistant from its center. Its spatial decay (shown in 
FigCH(b)) exhibits a power law decay with an exponent 
-1, consistent with earlier calculations [2^] in the free elec- 
tron limit. In the case of Cu, we used atom centered s— 
and c? — like orbitals and an s— like orbital centered at 
octahedral site O = (1, 1, l)f . Among atom-centered or- 
bitals, a d— like orbital (FiglSjc)) is quite localized and is 
occupied with 1.86e, while the s— like orbital is occupied 
with 0.69e. The orbital centered in the octahedral site 
(shown in Fig|3](d)) is occupied with l.Ole providing a 
6— centered bond, and exhibits slight mixing with atomic 
d— orbitals. Our results clearly reveal stronger directional 
bonding in Al than in Cu, which was also concluded from 
charge density analysis in understanding contrasting me- 
chanical behavior of Al and Cu 29] . 

Finally, we demonstrate application of the present 
scheme to localized orbitals in C2H6 and B2H6 molecules. 
In the former, we have used bond-centered orbitals and 
in the latter, we use bond-centered orbitals for the BH2 
radicals and s— like orbitals centered on H atoms that 
bridge the two BH2 radicals. C-C bonding orbital 
(FigBJa)) has a variance of 0.83 A 2 , which reduces to 
0.76 A 2 after joint diagonalization to obtain an MLWF. 
Orbitals corresponding to heteropolar C-H a bonds have 
their centroids closer to the H atoms. The orbital cen- 
tered on a bridging H atom of B2H6 (shown in Fig^Jb)) 
reveals the well-known three-centered bond, and has a 
variance of 0.96 A 2 each. 

In conclusion, we have presented a simple method for 
construction of well localized orbitals, that is applicable 
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